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We study the collective excitations in a dissipative charged fluid at zero chemical potential when an 
external magnetic held is present. While in the absence of magnetic Held, four collective excitations 
appear in the fluid, we Hnd Hve hydrodynamic modes in presence of magnetic Held. This implies 
that the magnetic Held splits the degeneracy between the transverse shear modes. Using linear 
response theory, we then compute the retarded response functions. In particular, it turns out 
that the correlation between charge and the energy Huctuations will no longer vanish, even at zero 
chemical potential. By use of the response functions, we also derive the relevant Kubo formulas for 
the transport coefHcients. 


I. INTRODUCTION 

Study of hydrodynamic limit in relativistic and non- 
relativistic systems has been an important issue since 
many years ago (See [I| and references therein.). In 
particular, study of transport phenomena in an external 
field, has attracted considerable attention recently. As a 
remarkable progress, it has been shown that in systems 
with chiral anomalies, the second law of thermodynam¬ 
ics necessitates the presence of parity odd terms in the 
constitutive relation of hydrodynamic current [H . 

To investigate the dissipative character of fluids, one 
way is to study the collective excitations around equi¬ 
librium in the system. In this paper we will consider a 
charged fluid at zero chemical potential and compute the 
dispersion relation of the collective excitations in pres¬ 
ence of an external magnetic field. We show that the 
presence of magnetic field splits the degeneracy between 
the transverse shear modes. While in the absence of 
magnetic field there are four independent hydrodynamic 
modes, we find five long wave modes here. Of these 
five, two are longitudinal sound waves and three others 
are purely dissipative shear modes. Moreover, compared 
to zero magnetic field limit, there appears new dissipa¬ 
tive contributions in all hydrodynamic modes. In other 
words, the presence of magnetic field intensifies the dis¬ 
sipation processes in a charged fluid. We will indicate it 
is actually the Lorentz force exerted from the magnetic 
field on the induced electric currents in the system that 
makes the fluid modes more quickly dissipate. 

As known, the microscopic information of fluids are 
reflected by the transport coefficients. These coefficients 
are in general related to equilibrium response functions. 
In order to determine how one could evaluate the trans¬ 
port coefficients in a magnetic field, we first compute the 
retarded response functions in presence of a magnetic 
field. Then, from S'0(3) covariant response functions, we 
derive relevant Kubo formulas for conductivity, shear and 
bulk viscosity as well. 

Since the presence of magnetic field breaks time re¬ 
versal symmetry in the system, the retarded functions 
G^j are no loMer symmetric under the exchange of 
”a” and ”b” Q. Specifically, the correlation func¬ 


tion between charge and momentum density fluctuations, 
namely G^.„, which was already vanishing, turns out to 
be non-zero and more importantly, anti-symmetric. We 
physically explain why the magnetic field forces the fluc¬ 
tuations in charge density to be correlated with those in 
momentum density. We also indicate, that the G^.„ is 
anti-symmetric, is the consequence of this fact that in 
presence of magnetic field, the symmetry under time re¬ 
versal holds if the sign of magnetic field is simultaneously 
reversed. 


II. RELATIVISTIC FLUID DYNAMICS 

In grand canonical ensemble at ^ yf 0, the equilibrium 
state is characterized by the density operator: 

z 

where Z = is the partition function. In the 

expression above, is the momentum operator, and N 
is the conserved charge operator. This means that in a 
general Lorentz frame, we need to know five parameters 
to specify the equilibrium state of a system: a tempera¬ 
ture T, a chemical potential ^ and three components of 
four-velocity = —1). 

Hydrodynamics describes the dynamics of a thermal 
system slightly deviated from its thermal equilibrium, or 
equivalently a locally equilibrated system, evolving to¬ 
wards its global equilibrium. To describe the hydrody¬ 
namic states, it is sensible to choose the slowly varying 
local functions T{x), p{x), u^{x) as the hydrodynamic 
variables. From the hydrodynamic equations which are 
simply the conservation laws of the conserved global cur¬ 
rents, one can find the solution for the hydrodynamic 
variables, just as a boundary value problem. 

In a general relativistic system with a C/(l) global sym¬ 
metry, the global currents are the energy momentum ten¬ 
sor and a current J^. So in presence of an external 
background field the hydrodynamic equations take the 
form 
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where the term in the right hand side of the first 

equation is the work which Lorentz force performs on the 
flow. Obviously, the number of unknown components (10 
components of and 4 components of J^) are more 
than the number of equations (4+1) [I^. However, the 
idea of local thermal equilibrium is that one can express 
and as functions of exactly five slowly varying 
hydrodynamic variables. 

Near the local equilibrium also, and may be 
expressed in derivative expansion of hydrodynamic fields. 
Up to first order in derivatives, the constitutive relations 
are given by: 

= [e + p)u^^u'' (3) 

= nu'^ + i/'^. (4) 


In the expressions above, the dissipative parts, namely 
and are not uniquely determined [l^. In the so- 
called Landau-Lifshitz frame Q we have: 


+ dpu^) - ( C - -77 ) P^"a.w(5) 


= -aTP>^’'d, 


(f) 




( 6 ) 


A. State of Equilibrium in presence of an External 
Magnetic Field 

Each equilibrium state of the system is specified by a 
static solution of the equations ([2]) . It is straightforward 
to see that the set {u^ = (1, 0), p = const., T = const.} 
not only characterizes state of equilibrium in the absence 
of external field, it can also be a thermodynamic solution 
when a constant background magnetic field is present. It 
is physically reasonable: magnetic field does not affect 
the static chargesfTsj. 

The case what we are interested in is exactly the above- 
mentioned equilibrium state but at /i = 0. Let us remind 
that a state in grand canonical ensemble at p = 0 has 
no charge, i.e. h = o[T^. In relativistic hydrodynam¬ 
ics, this means that there is an equal number of parti¬ 
cles and anti particles in thermal equilibrium. Hydrody¬ 
namic fluctuations in this state however would be a sen¬ 
sible thing, namely one can compnte Sp and 6n around 
equilibrium!^ . 

In summary, our goal is to compute the correlation be¬ 
tween hydrodynamic fields around the equilibrium state 

= (I, 0), T = const., p = 0, E = 0, B = const.. 


with = u^u’^+p^'^. Here, p, ( and cr are non-negative 
transport coefficients: shear viscosity, bulk viscosity and 
conductivity respectively. We may also define the elec¬ 
tric field in the rest frame of the fluid: . For 

consistency with derivative counting in hydrodynamic, 
we take the strength of the background gauge field 
of the same order of the temperature and the chemical 
potential; so ~ 0{d°) and Ffj_i, ~ 0{d) 0. This as¬ 
sumption constrains the value of magnetic field in our 
problem; it turns out that for hydrodynamic perturba¬ 
tion with momentum k, the magnetic field must be of 
order of H ^ yjp/ak. Such magnetic field could affect 
the velocity of sound in magnetized fluid. Since our sys¬ 
tem is assumed to be non-magnetized here, the effect 
of magnetic field appears at least from the first order 
in derivative expansion. We expect the magnetic field 
makes contributions of order of ^ p\f? to dissipative 
processes. 


III. HYDRODYNAMIC FLUCTUATIONS AT 
ZERO CHEMICAL POTENTIAL 


In order to study the hydrodynamic fluctuations we 
have to determine the state of equilibrium and then lin¬ 
earize the equations around that. Before proceeding fur¬ 
ther, let us denote that in what follows, rather than 
choosing T, p and m® as the hydrodynamic variables, we 
will choose their conjugate variables, the energy density 
e, momentum density tt® and charge density n. The im¬ 
portance of this choice is that the latter have microscopic 
definitions given by operators T°®(a;) and J°{x). 


B. Hydrodynamic Modes in Presence of Magnetic 

Field 


Having specified the state of equilibrium, we have to 
linearize equations around that. Let us remind that in 
the absence of magnetic field at /2 = 0, the fluctuations 
in charge density decouple from the fluctuations of mo¬ 
mentum and energy density (see ([5]). However when B 
is present in equilibrinm, even at p = 0, the fluctuating 
fields n(t, k), 7ri(t, k) and de(t,k) are coupled to each 
other as the following: 


dtn + Dn + —kjF^^iVk = 0 
w 

dtiT^ + ikjvlSe - = WF^^kmU + 

w 


dtSe + ikjTT^ = 0 


(7) 

with A4y = - kikj) - 'jshkj, w = e + p, "fr, = 

pIw, Js = {^P + C) /u!, D = a lx and x = {dnldp)^=Q 
The velocity of sound is defined through v1 = dpjde [^ . 

Equations (l7|) can also be written in a concise form. To 
proceed let us define the supefield (^a(Lk) = (n, tt®, e). 
In terms of 4>a^ the linearized equations can be written 
as 


dt(f)a{t, k) + Mab(k) k) = 0 (8) 

where the matrix M is 

/ k^P 0 \ 

Mab = -iDF^^km -M] - f PjmP™ 

\ 0 ikj 0 / 


( 9 ) 
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Now, hydrodynamic modes may be simply found via solv¬ 
ing the following equation: 

det (^-iujSab + Mab{k)j = 0 ( 10 ) 


The magnetic field itself, exerts a Lorentz force on the 
current or equivalently on the unit volume of the fluctu¬ 
ating element of the fluid as being 


Doing so, we obtain five collective modes around the equi¬ 
librium (See table (|T|). There are some interesting points 
regarding with these modes which we point them out 
in what follows. First, when B = 0, the modes corre¬ 
sponding to the momentum fluctuation in the directions 
transverse to k are degenerate, because the only preferred 
direction in the system is identified by the wave vector k. 
However, as expressed in table ll|^ , in presence of mag¬ 
netic field, the degenerate modes, namely UJ 4 and ws, are 
split into two dissipative modes. Moreover, when B || k, 
the degeneracy is restored again 

W4(k) = W5(k) = -i (y^k^ + |b2) . (11) 

The point which has to be physically explained is how the 
magnetic field participates in dissipation processes. Let 
us briefly explain how it happens for the case of shear 
modes. Consider a spacial frame wherein k = (fe,0, 0) 
and B = {Bx,By,0). In this frame, the shear mode W 5 
reflects the decay rate of the momentum fluctuations in 
the direction y. So we consider the momentum fluctu¬ 
ation as being tt^. According to ([5]), a momentum fluc¬ 
tuation in an external magnetic field, induces an electric 
current: 

Jz = -^B^TTy. ( 12 ) 

w 


dF= 3 xB= ( 0 , -^B^TTy, 0 ) . (13) 


No matter how the profile of Wy is, the above force is 
always directed opposite to iTy. Therefore, the Lorentz 
force makes the momentum of fluid element dissipates 
and the rate of dissipation is exactly the same as mag¬ 
netic contribution in With similar arguments, one 
can physically interpret the contribution of magnetic field 
to dissipative parts of the other hydrodynamic modes. 

As the second point, let us give a comment on the 
velocity of sound modes given in table HI As it can be 
clearly seen, the velocity of sound in a charged fluid is 
not affected by the presence of a background magnetic 
field. This is not surprising since the fluid is supposed to 
be non-magnetized here. If the fluid were magnetized, 
the pressure would not be isotropic even in equilibrium 
and therefore the sound velocity would be dependent 
on the direction of magnetic field. Investigating such 
effects would be of serious importance when studying 
magnetohydrodynamics in a charged system in presence 
of a strong magnetic field [3, @. 


Hydrodynamic modes in presence of magnetic field 

Hydrodynamic modes at B = 0 

a;i,2(k) = ±Vsk - f (k^y, + f B" sin^ 6) 

? 

u.i,2(k) = ±Vak - fk^ys 

0.3 (k) = , 

0.4 (k) = -iy,,k^ 

0.5 (k) = -iy,,k^ 

0.3,4(k) = -1 (^k2(D + y,) -p 1 B^ ± ^(^iD - y,) - f B^)2 + M^B^ k^ sin^ ej 

0.5 (k) = -i (k^y^ + :|B^ cos^ 6) 


TABLE I: Comparison between hydrodynamic modes ”in presence” and 
”in absence” of external magnetic field at /r = 0. In the expressions 
above, 0 is the angle between momentum vector k and the magnetic 
field B, i.e. cos 6 = fc.B. 


C. Response functions might be computed in the framework of linear response 

theory via using the following formula 


Using the methods of linear response theory, the re¬ 
sponse functions have been extensively discussed in non- 
relativistic hydrodynamics in Q. Recently, Kovtun ap¬ 
plied the methods of linear response theory to the case 
of relativistic hydrodynamics. Specifically, in Q, it has 
been shown that the retarded function 


k) = - (^dac + iooiK ^)ac) Xcb- (15) 

In the expression above, Kab = —iojSab + Afaf,(k) and 
Xab is the static thermodynamic susceptibility which is 
defined by 


Gfb(i-t',x-x') =-f6»(t-t') ([(^a(i,x),(^b(t',x')]^ (14) 


Xab — 


d(j)a 

dXb 


(16) 
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Here Xa is the external source term which couples to 
field (fia through the perturbation Hamiltonian. In or¬ 
der to identify sources, we have to find the linear change 
in Hamiltonian, namely SH{t,x.), for infinitesimal slowly 
varying disturbances (5r(t,x), 5/r(t,x) and ?;*(t,x). To 
proceed, we may formally write 


P(x){u^(x)P>^(x)+ iJ.{x)N(x)) _ -P{H+ SH(x)) 

^ z z ' 

(17) 

Simplifying equation (ED yields: 

SH{x) = - J tfxXa{x)^a{x) = 

— J (fix -I- Sfi{x)J^{x) + Vi{x)T^'‘(x)^ 

(18) 


Therefore we identify Xa = [6fi,Vi, as the 
sources corresponding to the hydrodynamic fields t/fo = 
{5n,'Ki,5e). As a result, the susceptibility matrix (IT6l) is 


(Tm 


Xab — 


ti/T 

0 w 




IJ./T 



/x 0 0 \ 

0 w 0 

Vo 0 f; 

(19) 


Similar to the case where the magnetic field is absent Q , 
the susceptibility matrix is diagonal here. This is due to 
this reason that chemical potential in our system is zero 
in equilibrium^. 

Now we substitute (m into formula (ESI to compute 
the retarded functions in presence of a constant external 
magnetic field. Results have been reported in table (HD, 
organized in terms of 50(3) group representations. In 
the scalar sector there are two non-vanishing response 
function. In oreder to represent the vector type response 
functions, we first have to find a vector basis. Using k 
and B, one can consider three independent vectors 

Vi=k, V2 = cos6 »B, V^ = Bxk (20) 

where cos0 = B.k. Taking these three vectors as the 
basis, we have covariantly expressed the vector retarded 
functions, namely and in the second row of 

table (HD- 

The basis of symmetric tensors may be constructed by 
using Sij plus all symmetric bi-vectors made out of basis 
vectors given above. It simply turns out that there are 
six independent symmetric tensors as being 

k,kj, m^n- ( 21 ) 

however only four of the them have participated in the 
expression of G^.^.. 


SO{3) Type 

Retarded Green’s Functions 

scalar 

Gf„(c.,k) = 0 ,Gf,(c.,k)= Gl(a;,k)= — 

vector 

1 

G«„(c^,k) = -G«^,(cc,k) = ^ (B X k), 

G«,(a.,k) = G«,(a;,k) = ^ {w{uj + ki + ia{B.k)Bi) 

Tensor 



+ ( + A345 ) + ( A1355 A345 ) 


TABLE II: Returded functions organized in different sectors of SO{3) 
group. In above expressions Ac = cos^ 6 -|- 7 ,,k^ and by the Aij... 
in denominators, we mean (cj — a;i(k))(a; — a;j(k)) .... 


In the absence of magnetic field, retarded response 
functions are symmetric under the exchange of a and 
b fields 

G«(cc,k) = G«(cc,k). (22) 

The origin of the above symmetry is that in a sys¬ 
tem which is microscopically time-reversal invariant, the 
space-time translational invariance implies thatQ 

Gfb(w, k) = 7?a??&Gf„(w, -k) (23) 

with rja being the time-reversal eigenvalue of the hydro- 
dynamic field (j)a ■ It is straightforward to show that when 


B = 0, das]) leads to ([22l).J^ 

When B is non-zero, and may participate in the 
structure of vector response functions too. On the other 
hand, the relatio n (1231) has to be generalized in presence 
of magnetic field 

G«(w,k,B) =ry„r?,Gf,(u;,-k,-B). (24) 

Denoting these considerations, it turns out that only 
is allowed to participate in the structure of G^. (See 
table (|ni))[2|. 

The situation is somewhat different in the case of G^.„. 
Consider a momentum density fluctuation appeared 
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at t = 0 around the origin. In a magnetic field being 
along the +y direction, the Lorentz force exerted on the 
fluctuating element, induces an electric current in the 
+z direction. Consequently, there appear fluctuations in 
density n at the points near the origin on the z axis. This 
means that in presence of B = (0, By,0), tt^ is correlated 
with n with the wave vector k=(0,0, ^ 2 ). As a result, 
the vector response function can only be made out 
of = B X k. The point which distinguishes from 
the other retarded response functions is that the former is 
not symmetric under the exchange of tt^ and n. It can be 
simply justified once one applies (l24l) to There are 

three minus signs which have to be taken into account: 
one from 77 ^. = — 1 , one from k —>■ —k and the last one 
from B -> -B. So 

G« Jc.,k,B) = -G«,^(c.,-k,-B) = 

(25) 

where in the last equality we have used the fact that 
GS^{uj, k, B) = scalar(B x k). 

It should be also noted that by taking the limit B —0, 
one can simply reproduce the results of Q. 

D. Kubo Formulas in Presence of a Magnetic Field 

Kubo formulas for transport coefficients of a charged 
fluid are well known in the literature^. However one 


IV. OUTLOOK 

Our results reported in this note might be encountered 
with a very important follow-up question. The question 
is how the hydrodynamic modes would be affected if we 
took into account the chiral effects in the fluid. Very 
recently, Yamamoto has shown that in a system with 
gravitational anomaly, there exists a new type of collec¬ 
tive mode. As shown in m , the so-called Chiral Alfven 
wave is a transverse hydrodynamic mode propagating in 


may be interested in computing transport coefficients 
for the charged fluid when existing an external magnetic 
field. Clearly, the response functions are changed com¬ 
pared to B = 0 case. Therefore the Kubo formulas are 
expected to be modified too. In table Ell we have listed 
the relevant Kubo formulas in presence of a constant 
magnetic field. 

It is straightforward to check that the B 0 limit 
of formulas given in table m coincide with the Kubo 
formulas for a charged fluid at zero chemical potential in 
Q. Let us emphasis that our Kubo formulas have been 
expressed in terms of response functions of hydrodynamic 
densities n, tt^ and e. However, one can also write them 
in terms of response functions of spacial currents. To 
proceed, one may first demand the response functions 
satisfy 


k^Gj^T,A^.k)= 0 (26) 

Gr^^jA^.k) = 77.„F“^Gj,j„(w.k). 


Then, by use of (1^ and (IMl) . it would be a simple cal¬ 
culation to exchange the response functions of densities 
with those of spacial currents in Kubo formulas. 


an incompressible fluid, with a velocity proportional to 
the anomaly coefficient. On the other hand, analogous to 
what we have done in this paper, one can study the dis¬ 
sipative effect in a chiral fluid. Specifically, it would be 
interesting to investigate whether the Alfven mode splits 
into two or more chiral modes when taking into account 
the dissipative effects or not. 

In another direction one may study the hydrodynamic 
modes and response functions in a charged Lifshitz fluid 
[Hill- We leave more investigation on the issue to our 


Retarded Green’s Functions 

(T = lim -^ 

k-»o k 2 + „ 

(“>k) 

7^ = 2 ^ lim 


(Bxk)’^ 2a^ 

X 

7 . = lim (Pif IniG« ..) - ^^2 

4(Bxk)^((B.k)2B2<T2-|-is2ij2) 2 2 

M T>s 


TABLE III: In above expressions we have used Vy = {B‘*‘a‘‘ — w^u^). 
My = (B'^cr^ -I- Vs = My [My — 2(B X k)'^w'^a'^. 


The projective tensors are defined as: 

= Sij — csc^ 0 kikj — csc^ OBiBj + cos 6 csc^ 9 k(iBj) 
= cot^ 9 kikj + csc^ 9BiBj — cos 9 csc^ 9 k^iBj) 
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future works. 
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